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Abstract 

Within the Local Potential Approximation to Wilson's, or Polchinski's, exact renormaliza- 
tion group, and for general spacetime dimension, we construct a function, c, of the coupling 
constants; it has the property that (for unitary theories) it decreases monotonically along 
flows, and is stationary only at fixed points — where it 'counts degrees of freedom', i.e. is 
extensive, counting one for each Gaussian scalar. Furthermore, by choosing restrictions to 
some sub-manifold of coupling constant space, we arrive at a very promising variational 
approximation method. 
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1. A C Function. 

Following Zamalodchikov's celebrated c-theorem[[l] for two dimensional quantum field 
theory, a number of groups have sought to generalise these ideas to higher dimensions ||. 
The motivation behind this, is not only to demonstrate irreversibility of renormalization 
group flows, in some precise sense and under certain general conditions, and thus prove 
that exotic flows such as limit cycles, chaos, etc. are missing in these cases, but perhaps 
more importantly to provide an explicit, and useful, geometric framework for the space of 
quantum field theories. In ref.|]T[|, Zamalodchikov established three important properties 
for his c function, 

1. There exists a function c(g) > of such a nature that 

| cs /)«(,) Acto)<o , (i) 

where the g % form an (infinite) set of dimensionless parameters, a.k.a. the coupling 
constants, and the beta functions are defined as (3 l = dg l /dt. The equality in (|IJ) is 
reached only at fixed points of the renormalization group g(t) = g*. 

2. c(g) is stationary at fixed points;0 i.e. (3 l (g) = for all i, implies dc/dg l = 0. 

3. The value of c(g) at the fixed point g* is the same as the corresponding (Virasoro 
algebra) central chargeH. (This property thus only makes sense in two dimensions.) 



Within the Local Potential Approximation0-[||(LPA) to Wilson's|3| or Polchinksi'sfllO 
exact renormalization group (which we describe below), we display a c-function which has 
the first two properties in any dimension D. This follows by virtue of the property that 
the /3-functions are the 'gradient flows' of c with respect to a positive definite metric. This 
property in turn follows straightforwardly from an observation by ZumbachQ. At fixed 
points, our c function is extensive, viz. additive in mutually non-interacting degrees of free- 
dom, as is also true of the Virasoro central charge (c/. property 3). Furthermore, we may 
normalise so that our c counts one for each Gaussian scalar and zero for each infinitely 
massive scalar (corresponding to a High Temperature fixed point), as does the central 
charge in two dimensions. It is probably not possible within the Local Potential Approxi- 
mation, to establish a more concrete link to Zamalodchikov's c. Of course, it would be very 
interesting to understand if these observations generalise to higher orders in the derivative 



1 We will not here need to restrict to critical fixed points [] 



expansion[ 11 ] [9J (which likely would allow a direct comparison with Zamalodchikov's c), 
or indeed generalise to an exact expression along the present lines. 

In the second part of this letter, we compute exactly within the LPA, some simple 
illustrative examples. In the third and final part, we point out that this geometrical 
structure for the LPA leads naturally to a variational approximation scheme, by restricting 
the flow to some finite dimensional submanifold of coupling constant space, of our choosing. 
Approximate continuum limits (i.e. fixed point behaviour), follow from stationarizing with 
respect to this finite parameter set. We briefly investigate the efficacy of the method, 
considering, in particular, the simplest polynomial approximation to the non-perturbative 
fixed point potential for a single scalar field in three dimensions. The resulting form, and 
critical exponents, lie exceedingly close to the exactH answers. 

We work with N Bose fields ip a in D Euclidean dimensions. We have checked that 
the constructions can be generalised to Fermi fields, but for simplicity we omit this. In a 



condensed notation|L2| [^3[, Polchinski's form of the Wilson renormalization group is given 

byraran 

dSA = 16SA dA uv SS A 1 dA uv d 2 S A 
dA 2 5ip a OA ' 5ip a 2 Y dA ' 5ip a 5y a ' [ ' 

Here S\[<p\ is the interaction part of a Wilsonian effective action S A = ^ a .A^.(^ a + 5A, 
and Ajjv = Cuv(q 2 / A- 2 )/q 2 is a massless propagator whose momentum q is cutoff by an 
effective ultra-violet cutoff A. We require that the cutoff function Cuv is analytic at 
q = and that Cuv(0) = 1 (so that physics is unchanged at scales much less than A), 
and Cuv ~ > as q — ► oo sufficiently rapidly that all momentum integrals in (|j) are 
well-defined (i.e. well regulated). Wilson's flow equation|3|] is identical to (0), after the 
transformation [1T5|1 [01 



<P 



a/ Cuv fi and identification Tt = — Sa • (3) 



In the Local Potential Approximation, we restrict the interactions to that of a general 
potential: Sa = fd xV((p, A), and discard from the right hand side of (§) all terms that 
correspond to higher derivative corrections (on expansion of Cuv where necessary). Thus, 

d a f dV \ 2 . D _ 3 d 2 V 



■%A) = - — - 7 A 



dA A 3 \d(p a J d(p a d(p a 

f d D q 

where a = -C' uv (0) and 7 = - J — —-^ C' uv (q 2 ) 



2 within LPA! 
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We assume that the dimensionless coefficients a and 7 are positive; this is the case, for 
example, if Cuv is strictly monotonically decreasing. To express all quantities in a form 
appropriate for the continuum limit (i.e. approach to fixed points), we change to dimen- 
sionlessi variables ip t— > iph. D l 2 ~ x */7, V 1— > VK D ^/a, and t = ln(/x/A) (where is some 
arbitrary physical mass scale), and by thus also absorbing the coefficients a and 7, obtain 
the LPA of Polchinski's equation 

d 1 dV d 2 V fdV\ 2 

^V(^t) + -(D-2)^ a ^-DV = -^—~(^-) , (4) 

at 2 dip a Oipadipa \d<p a ) 

in a manifestly 'scheme' independent form[|7|]. Since the LPA corresponds to discarding all 
momentum dependent terms from (Q), it is immediate to realise from (|3]) that the above 
is also the LPA for Wilson's equation. 

If in place of V, we introduce a Gibbsian-like measure p(<p,t) = exp{— V((p,t)}, and 
G(<p) = exp { — \(D — 2)<£> 2 }, and define the functional 



2 + ^p 2 (l-21np)j. , (5) 



then (^) may be rewritten manifestly as a 'gradient flow' 0: 

N dp bT 
a G m=-J-p ■ (6) 

In (||) and (^), a > is a normalisation factor for the ip measure, to be determined later. 
[Alternatively, changing variables ip 1— > <p/a, the factors a N can be absorbed and appear 
instead as factors of a 2 in front of the (d/dip) 2 terms in (f|) and (H).] Substituting the fixed 
point equation SJ-'/dp = back into (||), we see that fixed points p(cp, t) = p*(<p) satisfy, 

W = ja N jd N ^Gpl . (7) 

Now let g l (t) be a complete set of parameters (a.k.a. coupling constants) for V. Away 
from fixed points, these will be infinite in number. Asymptotically approaching a particular 
fixed point, it is possible to identify a finite number of parameters (namely the relevant 
and marginal couplings) and express all the rest (the 'irrelevant' parameters) in terms of 



3 



in general using A raised to the power of their scaling dimensions, however in LPA there is no 



anomalous dimension for (p. This fact also allowed our scaling-dimensionless choice for Cuv [ 15 1 . 
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them. (This 'reduction of parameters' corresponds to the continuum limit defined around 
this fixed point, see later and e.g. |P2]] ||16|1 ; universality arises because perturbations in the 



irrelevant parameters about these values, exponentially decay away as t — > oo.) Following 
ref. |]]], we define the operators 

<$> i {g) = d i V{ip,g) , (8) 

where we have written di = d/dg l . We define the metric 

g iJ (g)=a N Jd N ^Gp 2 ^ J . (9) 

It is positive definite providing V is real, which we assume (corresponding to a unitary 
theory in Minkowski space). Multiplying (|6|) by dip = —p&i, we thus obtain 

GyP = -diF(g) . (10) 

We define the c(g) function through 

F = DA c /4 , (11) 

where A > 1 is a normalisation factor to be determined. Redefining the metric by a 
positive factor, = (.Fin A) Qij, we obtain from (p7Q|) : 

d i c(g) = -g ij P(g) . (12) 

From (|T2| ) and the positive-definiteness of Qij, we see that ([!]), and the properties 1 and 
2 hold. Now suppose that the fields form two mutually non-interacting sets. Let us write 
tp a = (fa when the field belongs to the first set, and tp a = when it belongs to the 
second set. Similarly, the couplings g split into two sets g^ and g^ 2 \ The potential V(ip, t) 
can be written V(<p,t) = V^{(p^\t) + V {2 \ip^ 2 \ t). Thus p factorizes: p = p (1) p (2) , and 
from (0), at fixed points we have J-"[p*\ = j}J-[p* ]J-[p* ]. Therefore from (|Tl|) we have, 
at fixed points, that our c is extensive: 

c(g*) = c(gP) + c(gl 2) ) . (13) 

At the Gaussian fixed point V* = 0, the Virasoro central charge counts one degree 
of freedom per scalar, i.e. is here equal to N At the High Temperature fixed point 
V* = V^ T = hip 2 — N/D, the Virasoro central charge vanishes, because V^ T corresponds 
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to a non-critical fixed point, i.e. infinitely massiveu fields and thus no propagating degrees 
of freedom. This interpretation of V^ T follows directly from the Legendre transform 
relation between the Polchinski and Legendre effective potentials |13|] [|12], or by the explicit 
flow derived below (see also ref. ||17|1). 

The normalisation factors A and a may be uniquely determined by requiring that 
our c agree with this counting in any dimension D. Substituting c = in (|Tl| ) and 
p* = exp — V^ T in (0), we obtain 



-2/D D + 2 
• = e V— ' (14) 

For the Gaussian fixed point we substitute c = iV in (^TJ) and p* = 1 in ([7|), and thus 



-2/D [D + 2 , . 

A = e ^j— . (15) 

Note that A > 1, as required, for all D > 2. The precise values for a and A may be 
expected to change beyond the LPA. 

Perturbing the couplings to first order about a fixed point: g l (t) = g\ + ev l e , we 
obtain from (|T2|) the eigenvalue equation 

didjc{g*) v J = -\Gij(g*) v J , 

determining the eigenoperators $ = v l &i and their scaling dimensions D — X. Equivalently, 
from (0), 

didjJ^ig^) v J = -\Gij(g*)v J . (16) 

2. Examples. 

Consider the simple example of the Gaussian fixed point perturbed by the mass oper- 
ator, for a single scalar field. Thus we set V((p, t) = ^a(t)(p 2 +£(t). The (3 functions for 
a and £ follow easily from (|4]): 

-a = 2a(l-a) 

dt (17) 
m 8 = D£ + a . 



in units of A 
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The general solutions are 

rr(+\ — 

-2t 



1 + r e 

S(t)=S e ~ ~ du ^ttt , (186) 



2 J i + ure 2t 

where r and So are integration constants. Considering t — > — oo, one sees that this solu- 
tion indeed emanates from the Gaussian fixed point, while a change of variables back to 
dimensionful (i.e. physical) variables shows that r ~ |U 2 /m 2 , where m is the mass of the 
still-Gaussian scalar field, while So[± D is an added vacuum energy term. We normalised 
the special solution in (|18£| ) so that with £ = 0, (|18| ) tends to V^ T as t — > oo. 
Since p is Gaussian, @ is readily determined; 

_ a -25 f 6Da + 4a 2 + D 2 - 2D) I n , . 

T=-n < 2D£ H L/ . (19) 

Combining this, ([14]), ([15]) and (|TT1), yields c(t) and one verifies that when So = 0, c(£) 
flows from 1 to as t runs from — oo to +oo. Note that if So ^ 0, then c(t) — > =Foo as 
t — > oo, depending on the sign of £q- This seems in contradiction with the idea that c 
counts degrees of freedom, since evidently the vacuum energy should not figure per se in 
this counting. However, one must recall from (|TB|), that c is extensive (i.e. 'counts') only 
at fixed points, and with So ^ the system never reaches another fixed point as t — > oo. 

From the Gaussian form of G, we recognize that at both the Gaussian and High Temper- 
ature fixed points, the metric (H) (and thus Qij) is diagonalized by choosing the operators 
$i to be products of Hermite polynomials H n in the <p a . Since these also turn out to 
diagonalize didj\F(g*), they are the eigenperturbations, and the corresponding eigenvalues 
follow straightforwardly. Choosing the Gaussian fixed point = 1 for example, and again 
specializing to the case of one scalar field for simplicity, we thus take $ n = H n (f y/D - 2), 
n = 0, 1, • • •. The metric has non-zero components ||18|| Q nn = 2 n+1 n\a^j D n _ 2 . From @ 
and fl5[), we obtain 



didjT 



<9$ <9<1> 

-D +a ld(p G 3 



P =i 13 "J " r d(p dip 

Thus, using! H' n = 2nH n _i, we have that didjT is also diagonal, and we recover the 
expected Gaussian spectrum of eigenvalues A = —d n d n J-'/Q nn = D + i(2 — D)n. 



prime being differentiation with respect to the argument 
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We may also construct perturbatively, the continuum limit about the Gaussian fixed 
point, and c and Qij, to any desired order in the relevant and marginal couplings. Thus 
for example in D = 4 dimensions, the massless continuum limit is constructed by solving 
the j3 functions for g 2 , g 6 , g 8 , etc. , iteratively in terms of a power series in g (t) (in direct 
generalisation of the case for "A(t)" and "71 (£)" given in ref.||12||), reducing ( p!2|) to the one 
beta function /3 4 (g 4 ) = — d^c^) /Q^. 

3. Variational Approximations. 

The gradient flow form (|(J) of the LPA, suggests the possibility of approximating p by 
a variational ansatz, that is setting p = p((p; <7i, • • • , (?m), M < 00, where p is some finitely 
parametrized set of functions of our choosing. Interpreting the functional derivative in (Q) 
to include only variations in this restricted set, we arrive again at ([10]) [or (|T2|)1 where 
however, here and from now on, indices run only from 1 to M. Thus geometrically, we 
restrict the flows to the sub- manifold M. parametrized by gi, ■ ■ ■ , gu- 

Approximations to the fixed points of (|]) follow, by (|10D, from solutions to the vari- 
ational conditions diT = 0,i and correspond geometrically to those points g = g* on M. 
where T g M. (the tangent space to M. at g) is perpendicular to the exact! flows. From this, 
it is intuitively clear that good results can be expected generically, e.g. for the shape of p, 
if M. passes sufficiently 'close' to the true fixed points, and good results will be obtained 
for the scaling dimensions and shape, of any such eigenoperators that are almost parallel 
with Tg^Ai. We can expect decent approximations for the full flows, from say, one fixed 
point to another - under similar conditions. 

It is clear that it is possible to be unlucky and find in this way 'spurious fixed points' 
that do not well approximate the exact solutions, but these will not be relatively stable 
under changing or improving the ansatz manifold Ai, e.g. by enlarging its dimension M. 
In fact as M. improves, such spurious fixed points (if indeed there are any) will disappear 
entirely. Equally, it is clear that numerical results can (at least in principle) be obtained 
to any desired degree of accuracy. 

These characteristics are in marked contrast to the general situation for truncations (of 



the renormalization group) to a finite set of operators |T^]0], where spurious fixed points 



multiply at higher orders, and numerical results cease to converge. 



6 since the (restriction of the) metric Qij is still positive definite 

7 i.e. with M = 00, or from (13) 



To test the practicality of this variational method, we tried the following very simple 
even polynomial ansatze. The simplest two cases, V(ip,t) = £{t) and V(ip) = £{t) + 
^a(t)v? 2 , both correspond to exact reductions of the full flow equations (a.k.a. (3 functions) 
from (Q), and are given by fll~7T). Therefore it is immediate to realise that we obtain only 
—and exactly — the Gaussian and, when a is included, High Temperature fixed points. 
Similarly, we obtain the eigenvalue A = D for the unit operator $£ = dgV = 1 at both 
fixed points, and A = 2 (-2) for the second eigenvalue at the Gaussian (High Temperature) 
fixed point when a is included. This may be checked directly through ([19]) , (|l~6l), and (||). 

In fact it is straightforward to confirm that dgV = 1, A = D, is always a solution in any 
approximation that includes £ as a parameter. Equally, from our previous exact analysis 
of the Gaussian and High Temperature fixed points, it is clear that a general polynomial 
ansatz of order <p n , n > 0, with unconstrained coefficients, will find exactly these fixed 
points and exactly the first n eigenvalues. 

The next-simplest even polynomial is V((p, t) = £{t) + ^a(t)ip 2 + s(t)(p 4 . Since the 
case s = has already been analysed, we fix s > 0. We set D = 3. We take advantage 
of coordinate invariance to write s = (g 4 ) and a = g 4 g 2 — 1/4, and change variables 
V 9 = x l Vy ) so that T and its derivatives may be expressed in closed form in terms of the 

integrals I n (9 2 ) = j_ dx x n e , n = 0,2. In this way, on the basis of analytic 

(for large ±g 2 ) and numeric estimates, we establish that there is only one non-trivial 
solution of the variational equations diT = 0, i = £, 2, 4. It corresponds to s = .00772624, 
a = —.13488 and £ = .054794. In fig.l, we plot the resulting form for and compare it to 
the one exact non-trivial fixed point solutioni from (Q) (which itself is an approximation 
to the Ising model fixed point in three dimensions). Solving (^6|), we obtain, apart from 
A = 3, u = 1/X = .6347 from the positive eigenvalue, and u> = —A = .6093 from the 
remaining eigenvalue. These are 2% and 8% off the exact values from (^), namely v = .6496 
and iv = .6557, respectively [0. Overall, we find these results, for this simplest possible 
variational approximation to the non-trivial fixed point (and/or u), truly impressive. 

Since the fixed point behaviour of (|j) for a single field — or more generally — a single 
invariant, is straightforward to solve directly numerically, the true potential of this vari- 
ational method lies in the relative ease in which approximations for global flows may be 
solved, and approximate solutions of LPAs found for more than one invariant, which thus 
correspond to (possibly high dimensional ^tj) partial differential equations. 

Should the c function be extended beyond the LPA, we expect that the variational 
method will prove to be even more powerful. 

8 Solved by shooting as described in refs. |i~9|l [f7[|. 
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Fig.l. Plotted as p*, the simplest polynomial variational approximation to the non-trivial 
fixed point potential (dashed line) is compared to the exact solution to (full line). 
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